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Abstract. Dynamic Spectrum Access systems exploit temporarily available spectrum ('white 
spaces') and can spread transmissions over a number of non-contiguous sub-channels. Such 
■ methods are highly beneficial in terms of spectrum utilization. However, excessive fragmentation 

^v^j degrades performance and hence off-sets the benefits. Thus, there is a need to study these 

processes so as to determine how to ensure acceptable levels of fragmentation. Hence, we 
present experimental and analytical results derived from a mathematical model. We model a 
system operating at capacity serving requests for bandwidth by assigning a collection of gaps 
"^J>< ' (sub-channels) with no limitations on the fragment size. Our main theoretical result shows that 

even if fragments can be arbitrarily small, the system does not degrade with time. Namely, the 
average total number of fragments remains bounded. Within the very difficult class of dynamic 
fragmentation models (including models of storage fragmentation), this result appears to be 
the first of its kind. Extensive experimental results describe behavior, at times unexpected, of 
fragmentation under different algorithms. Our model also applies to dynamic linked-list storage 
allocation, and provides a novel analysis in that domain. We prove that, interestingly, the 50% 
rule of the classical (non-fragmented) allocation model carries over to our model. Overall, the 
^ ' paper provides insights into the potential behavior of practical fragmentation algorithms. 
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1. Introduction 



This paper focuses on dynamic resource allocation algorithms in Dynamic Spectrum Access 
Networks (also known as Cognitive Radio Networks) . A Cognitive Radio is a concept that was first 
defined by Mitola [2"Tll2"B] as a radio that can adapt its transmitter parameters to the environment 
in which it operates. Technically, it is based on the concept of Software Defined Radio [4] that 
can alter parameters such as frequency band, transmission power, and modulation scheme through 
changes in software. According to the Federal Communications Commission (FCC), a large portion 
of the assigned spectrum is used only sporadically [1 lj . Because of their adaptability and capability 
to utilize the wireless spectrum opportunistically, algorithms for Dynamic Spectrum Access are 
key enablers to efficient use of the spectrum. Hence, the potential of Cognitive Radio Networks 



has been recently identified by various policy p~2|fT3| , research [318] , and standardization pT0]l 19|l20j 
organizations. 

Under the basic model of Cognitive Radio Networks [2j, Secondary Users (SUs) can use white 
spaces (also known as spectrum holes) that are not used by the Primary Users but must avoid 
interfering with active Primary Users (e.g., Figure [T]). For example, the Primary and Secondary 
Users can be viewed as TV broadcasters and cellular operators using available TV bands p"9] . 
Under this model, one assumes that when a transmission of a Primary User takes place, it occupies 
a predefined band. An SU may identify spectrum holes (not used by Primary or other Secondary 
Users) and can allocate its bandwidth among a number of subchannels, occupying a number of 
holes (not necessarily contiguous). This can be realized, for example, by employing a variant of 
Orthogonal Frequency-Division Multiplexing (OFDM) that is capable of deactivating sub-carriers 
which have the potential to cause interference to other users [16,25,29 32] (such a non-contiguous 
OFDM scheme is shown in Figure [T]). 

The use of non-contiguous bandwidth blocks results in non-trivial behavior even for very simple 
scenarios. As an example, due to the dynamic use of the available spectrum holes and the arrivals 
and departures of SU bandwidth requests, SUs may need to transmit in a set of smaller and 
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Fig. 1. Non-contiguous OFDM with Primary and Secondary Users (SUs), where 
the Secondary Users use non-contiguous channels that do not overlap with the 
Primary Users' channels. 
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Fig. 2. An example of the admission and departure processes and the resulting 
spectrum fragmentation. 



smaller holes. This will lead to a highly fragmented spectrum whose maintenance may require 
complex algorithmic solutions. Although the practical (physical layer) aspects of OFDM-based 
Dynamic Spectrum Access have been extensively studied recently, the use of fragmented (i.e., 
non-contiguous) spectrum introduces several new problems |21il31[[33] that significantly differ 
from classical Medium Access Control (MAC) and fragmentation problems. 

In this paper, we study the most basic theoretical model in which the spectrum is shared by 
SUs only. Those users have to transmit and receive data, and accordingly need some bandwidth 
for given amounts of time. Hence, bandwidth requests of SUs arc characterized by a desired total 
bandwidth and the duration of a time interval over which it is needed. The data transmission can 
take place over a non-contiguous channel (i.e., a number of subchannels). Once a transmission 
terminates, some fragments (subchannels) are vacated, and therefore, gaps (spectrum holes/white 
spaces) develop randomly in both size and position. 

When allocating a channel to a new SU request, it is being fragmented (in the frequency domain) 
into available gaps until the full requested bandwidth is provided. This process repeats itself, until 
the next request fails to fit into the available fragments. Figure [2] demonstrates the process of 
a transmission termination followed by the fragmentation of waiting bandwidth requests (more 
details about this example can be found in Section [2]). We note that from a practical point of 
view, such a system can be viewed as a simplistic version of an OFDM-based access-point/base- 
station/spectrum-broker that tracks the available spectrum holes and allocates non-contiguous 
bandwidth blocks to SUs. 

The main goal of the paper is to investigate the phenomena of fragmentation induced by spec- 
trum allocation algorithms. For this purpose, we will ignore the particularities of techniques such 
as OFDM and make a couple of assumptions (for a complete system description, see Section [2]): 
(i) the system operates at capacity and there is always a waiting bandwidth request; and (ii) the 
fragment size is not bounded from below. Making the first assumption allows us to study the effect 
of fragmentation in the worst case. Clearly, if there are idle periods, when there are no waiting 
requests, only departures occur and the fragmentation level of the system decreases during these 
periods. Similarly, the latter assumption allows us to study the system performance when artificial 
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lower bounds on the fragment size arc not imposed. This differs from OFDM-based systems in 
which a subcarrier has a given minimal bandwidth. 



Main Results. At first glance, it seems that given assumption (ii) above, as time passes, the 
fragments used by bandwidth requests might well become progressively narrower and that the 
number of fragments grows without bound. Clearly, such an operation model cannot be supported 
by a realistic system (e.g., an OFDM-based system). However, our analytical and experimental 
results show that this is not the case. 

In particular, under our model, a portion of the spectrum (represented by the interval [0, 1]) is 
available for allocation to requests whose sizes are uniform on (0,a] (0 < a < 1). Once allocated, 
requests are satisfied (leave the system) after an exponential time. A stochastic model of channel 
fragmentation is investigated which includes the asymptotic behavior of the number G(t) of gaps, 
F(t) of fragments, and R(t) of requests being served (also referred to as allocated channels). A 
Markovian description of this system is clearly more complicated than the vector (G(t), F(t), R(t)), 
since for each of the R{i) channels, the location and size of each of its fragments should be included 
in such a description. 

Let (t n ) be the sequence of request departure times. Our main result shows that the total 
number of gaps and fragments at these epochs has bounded exponential moments, 



for some r\ > 0. It follows that the sum F(t n ) + G(t n ) is strongly concentrated near the origin, 
indicating that, with high probability, low fragmentation of the spectrum holds under high traffic 
intensity. In our analysis, several basic ingredients are used: relations between the number of 
fragments and gaps in the spectrum; a drift relation of Lyapunov type for the total number of 
gaps, fragments, and requests (Proposition!!]); a general inequality for Markov chains (Theorem|4]); 
and a stability result of [22] . 

In addition to the analytical results, and in order to gain insight into the performance of the 
system, we present the results of an extensive simulation study. We show that although for a 
given maximum request size (a), the number of fragments has a finite expected value, there is 
a linear relationship between 1/a and the expected number of fragments into which a request is 
divided. This indicates that when the requests are very small, they are also fragmented into a 
relatively large number of fragments. From a practical point of view, this implies that imposing 
a lower bound on the fragment size is important in order to reduce the complexity resulting from 
maintaining a channel heavily fragmented. 

We show that different bandwidth allocation algorithms have significantly different performance 
in terms of the fragmentation occurring in the system. In particular, an algorithm (referred to as 
LFS) that allocates gaps in decreasing order of their sizes reduces the fragmentation by almost an 
order of magnitude compared to other algorithms. Interestingly, we also show that the number 
of fragments is distributed according to a Normal law with a relatively low mean value. Finally, 
we observe an unexpected reappearance of the 50% rule [23], according to which, on average the 
number of gaps (unoccupied fragments) is half the number of requests being served (this rule 
differs from the original 50% rule by Knuth [23] which is briefly discussed below). 

Related Work. The areas of Dynamic Spectrum Access and Cognitive Radio Networks have been 
extensively studied (for a comprehensive review of previous work see [Tl[2l l24[[2l)] ). In particular, 
problems stemming from fragmented spectrum have been considered in [21,31,33 . In addition, 
several previous works focused on the technical aspects of channel fragmentation 15, 16, 18,25, 
29,30,32 . To the best of our knowledge, despite these recent efforts, the fragmentation problem 
considered in this paper has not been studied before. 

Yet, the spectrum allocation problems described here can be characterized in the context of 
classical dynamic storage allocation problems of computers, see Knuth [23]. In that context, the 
term "spectrum" refers to the storage unit, "bandwidth" refers to storage, and "channel" refers to 
a region in storage containing a file. In the original storage model, fragmentation refers only to the 
gaps of unoccupied storage interspersed with intervals of occupied storage: files are not fragmented. 




n>l 
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The 50% rule was derived for this model and asserted that the number of gaps was approximately 
half the number of files when exact fits of files into gaps were rare. Our notion of fragmentation 
applies also to the files themselves, so our discovery that the 50% rule continued to apply was a 
surprising one at first. Note that, in terms of the storage application, our model corresponds well 
to linked-list allocation of files - channels allocated to requests are sets of linked, disjoint segments 
of storage allocated to specific file fragments. Our results provide a novel analysis for such systems, 
and have implications for the garbage-collection/defragmentation process in linked-list systems. 

We note that studies of dynamic storage allocation have been around for some 40 years, and 
widely recognized as posing very challenging problems to both combinatorial and stochastic mod- 
cling and analysis. In particular, results of the type found in this paper, rigorous within stochastic 
models, seem to be quite new. 

For a system without fragmentation, the early results of Kipnis and Robert |22] concern a 
stochastic analysis of the number of requests being served in the system. Fragmentation is not an 
issue in their case since request allocations can be moved when needed in order to put available 
space all together in one block. A major result in [22] applying to our model asserts the existence 
and uniqueness of an invariant measure for the number of requests. Explicit formulas for our 
system are hard to come by, but those in |22j for the maximal departure rate in special cases have 
provided useful checks on our experiments. 

Organization of the Paper. The remainder of this paper is organized as follows. Section [2] 
describes the system, allocation algorithm, and the mathematical variables describing fragmen- 
tation. Section [3] presents simulation results that bring out the behavior of the system and the 
effects of fragmentation. Section Q] introduces some relations between the variables describing the 
fragmentation and a 50% limit law for the relation between the number of active channels and the 
number of gaps in the spectrum at departure times. Section [5] contains our main mathematical 
result, Theorem [2 which shows that the average value of the total number of fragments and gaps 
is bounded. Furthermore, in some cases, the existence of an equilibrium distribution is established 
in Theorem [3] Section [5] discusses algorithmic issues, such as changes in performance resulting 
from alternative algorithms for sequencing through the list of gaps when constructing a channel 
for a newly admitted request. Section [7] presents experiments suggesting that Normal approxima- 
tions for the total number of fragments and gaps hold. Section |8] discusses the results and future 
research directions. 

2. The Model 

System Model. As mentioned above, we consider requests for bandwidth queued up waiting 
to be served, each of them identified with the amount of bandwidth required and a time telling 
how long it wants a channel with that total bandwidth. Channels are allocated to bandwidth 
requests on a FCFS (first-comc-first-served) basis, subject to available bandwidth; the channels 
must remain fixed while active; they depart after varying delays, so gaps alternating with sequences 
of sub-channels (channel fragments) develop over time. 

For convenience, we normalize the spectrum to the interval [0, 1], so all bandwidth requests are 
numbers in [0,1]. There is a queue of waiting requests that never empties. Assuming that the 
spectrum initially has no channels in use, the allocation process begins by allocating consecutive 
channels i.e., consecutive subintervals of [0,1], to requests whose sizes are u±, . . . until a request 
whose size Ui, i > 1, is reached which exceeds available bandwidth, i.e., 

u\ + • ■ ■ + Ui-i < 1 < Ui H + iii 

All i — 1 of the channels now begin their independent delays. Subsequent state transitions 
take place at departure epochs when the delays of currently allocated channels expire. At such 
epochs, all fragments of the departing channel are released. Suppose all requests up to Uj have 
been allocated channels and a request Ui, i < j, departs, releasing its allocated channel. Then 
Uj+x, Uj+2, ■ ■ ■ are allocated their requested bandwidths until, once again, a request is encountered 
that asks for more bandwidth than is available. All channels then begin or continue their delays 
until the next departure. A standard rule for setting up a channel is to scan the spectrum, 
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from one end to the other, with gaps of available bandwidth allocated to fragments of the new 
requested bandwidth until enough has been allocated to satisfy the entire request. The last gap 
used in satisfying a request is normally only partially used; the partial allocation in the last gap 
is left-justified in that gap. We refer to this scanning rule as Linear Scan. 

An example is shown in Figure [2] Allocations for the first 8 user requests whose sizes are 
Mi, • • ■ , Us are shown, assuming that, at time 0, the spectrum is not in use. The requests with 
sizes u\ 1 U2, and U3 are the first to be allocated channels; M4 must wait for a departure, since the 
first 4 request sizes sum to more than 1. The variables U give the sequence of departure times of 
allocated requests. We see that the first occurrence of fragmentation takes place at the departure 
of it2 and the subsequent admission of 1*4; an initial fragment of U4 is placed in the gap left by 1*2 
and a final fragment is placed after 113. Note also that, even after 1*2 and 1/4 have departed, there 
is still not enough bandwidth for 145. After the additional departure of iti, both 1*5 and 1/6, but 
not M7, can be allocated bandwidth. 

In Scction[6j we shall also evaluate two alternatives to the Linear Scan rule: a Circular Scan, by 
which each linear scan starts where the previous one left off, and a Largest-First Scan intended to 
further reduce fragmentation. Although these algorithms make different scans of the gap sequence, 
they are all alike in their treatment of the last gap occupied: the last fragment is left justified in 
the last gap. This is a key assumption, and it is very likely to hold in practice. 

The allocation process described in this section contrasts with the classical model of dynamic 
storage allocation, in which each bandwidth request must be accommodated by a single, sufficiently 
large gap of available spectrum. However, the process does correspond closely to the linkcd-list 
model of dynamic storage allocation, in which an available-space list is maintained and files can 
be fragmented in accordance with this list. This paper also yields a novel stochastic analysis of 
such systems. Note that the assumption that there is always another request in queue models a 
system operating at capacity, where the departure rate, which is equal to the admission rate, is 
often called maximum throughput. 

Notation, State Space, and Probability Model. We denote by U the size of the request 
waiting to be allocated bandwidth, and by Ui for i > 1 the size of the ith request behind it. 
Except in the proof of Theorem [3] in Section [5j we omit the dependency in time of these variables 
for the sake of notation. As mentioned in Section O we denote by R(t) the number of channels 
(requests that are allocated bandwidth) at time t. F(t) and G(t) denote the number of fragments 
and gaps at time t. A Markovian description of this system is clearly more complicated than the 
vector (G(t),F(t), R(i)). Hence, we now define the state space of the fragmentation process. At 
a given state, we denote by r the number of channels (requests to which bandwidth has been 
allocated) and by Si (1 < i < r) the amount of bandwidth allocated to these requests. A state of 
the spectrum [0, 1] carries the information given by a sequence in which gaps alternate with sets 
of contiguous fragments: 

Definition 1. A state x in the state space S of the fragmentation process is denoted by 

x = (L lt ...,L r ;u) 

where u is the size of (amount of bandwidth required by) the request waiting at the head of the 
queue and Li, is the list of open subintervals of [0, 1] occupied by the fragments of the i-th channel. 
For x to be admissible, the open intervals in UiLi must be mutually disjoint, and, since the size of 
u exceeds the bandwidth available, u > 1 — Sj has to hold. 

Note that for a specific state, the size of the requests in the queue are denoted by it and 
Uj , j > 1 and the size of the channels (requests that are already allocated bandwidth) are denoted 
by s u 1 < i < rQ 

Let (X(t)) be the process living in the state space S. In the probability model used in this paper, 
the bandwidth requests (Ui) arc independent random variables which are uniformly distributed 



For simplicity of the presentation, we did not use the notation Si in the example provided in Figure [2] 
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Fig. 3. Average numbers of channels, of gaps, and of fragments per channel. 

on (0, a],0 < a < 1, and request residence times arc i.i.d. mcan-1 exponentials^ Under this 
probability model, (X(t)) is a Markov process on S. 

3. Experimental Results 

In this section we describe an experimental study of the model described above. This serves 
two related roles. First, it brings out characteristics of the fragmentation process that need to be 
borne in mind in implementations, particularly where these characteristics show conditions (e.g., 
parameter settings) that must be avoided, if a system with fragmentation is to operate efficiently. 
The second role is that of experimental mathematics, in which results indicate where behavior 
might well be formalized and rigorously proved as a contribution to mathematical foundations. In 
the latter role, this section leads up to the next two sections, which formalize the stability of the 
fragmentation process. 

The experiments were conducted with a discrete-event simulator written in C. In general, the 
tool is capable of simulating stochastic request arrival/departure processes. For this paper, how- 
ever, the arrival process was effectively inoperative, since the interest here is behavior while the 
system is operating at capacity (i.e., at maximum throughput). Admissions are made whenever 
the waiting request fits into available spectrum, and once made the waiting request is immediately 
replaced by another with independent samples from the required-bandwidth and residence-time 
distributions. The admission process continues in this way, effectively simulating a queue that 
never empties, until no further admissions can be made and one or more departures need to occur. 
The excellent accuracy of the tool was established in tests against exact queueing results, and 
exact results from [22] . The verification details are omitted and can be found in Appendix lAl 

Recall that bandwidth requests arc uniformly distributed on (0,a],0 < a < 1. Hence, the 
principal parameter of the experimental model will be a. The simulations of stationary behavior 
were most demanding, of course, for small a. For every a value, 0.01 < a < 1, 20 million 
departure events were simulated starting in an empty state, with data collected for the last 10 
million events. For 0.001 < a < 0.01, 100 million departure events were processed and data 
collection was performed during the last 50 million events. 

The results for the average number of channels, the average number of gaps, and the average 
number of fragments per channel arc shown in Figure [3] The curves are nearly linear in 1/a (the 
errors in the linear fits are within the thickness of the printed lines). In particular, the asymptotic 
average number of channels in the spectrum is 2 /a. When requests are large relative to the 



The assumption that request sizes are uniformly distributed is for convenience. The results hold for more 
general distributions. 
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Fig. 4. Average total number of fragments vs. 1/a: a quadratic fit. 




Fig. 5. Percentage of type-i fragments. 



spectrum (i.e., for a > 1/3), the behavior is not given by functions quite so simple. As such cases 
are of less practical interest, we omit the relevant data due to space constraints. 

The asymptotic linear growth of the average number of channels as a function of channel size 
is obvious, but the linearity of the other two measures is not so obvious. A closer look shows 
that the average number of gaps is almost exactly one half the average number of channels for 
even relatively small 1/a. This is an unexpected version of Knuth's 50% rule for dynamic storage 
allocation. We return to this behavior in the next section, where we prove a 50% limit law. The 
linear growth of the average number of fragments per channel may also be unexpected at first 
glance: the fragmentation of channels increases as the average channel size decreases. This linear 
growth implies the quadratic growth of the average total number of fragments plotted in Figure [4] 
(the accuracy of the fit is as before: the error is within the thickness of the printed lines). 

The analysis in the later sections will focus largely on tracking fragment types defined as follows: 
a fragment is of type-i, if it is adjacent to 0, 1, or 2 fragments. It can be seen in Figure [5] that for 
small a, more than 90% of the fragments are type-2 fragments. In addition, clearly, the number 
of type-0 and type-1 fragments is a function of the number of gaps. These observations and the 
results illustrated in Figures [3] and 0] indicate that, even for relatively small 1/a, the average total 
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Fig. 6. Average sizes of fragments and gaps. 



number of type-0 and type-1 fragments grows linearly in 1/a, but the average number of type-2 
fragments grows quadratically. 

Figure [5] compares the average gap and fragment sizes. As might be expected, for relatively 
small a, they are close to each other. The relation holds even for moderately large a, although 
for a rather close to 1, the difference amounts to about a factor of 2. With this property and the 
50% rule suggested by Figure El the linear growth in the number of fragments per channel (shown 
in Figure [3]) is easily explained for moderately small a in the following way. 

As mentioned above, for moderately small a, the number of channels is approximately 2/ a (i.e., 
the spectrum size divided by the average request size). Due to the 50% rule, the number of gaps is 
roughly 1/a. At any time, the total size of the gaps is at most a, since there is a request waiting 
for departures whose requested bandwidth exceeds that total size of gaps. Therefore, at most a 
available bandwidth is spread among 1/a gaps, giving an average gap size on the order of a 2 . The 
average fragment size is at most (and indeed very close to) the average gap size. The fragments 
must occupy at least 1 — a of the spectrum (since, as mentioned above, at most a fraction a of 
the spectrum is devoted to gaps). Thus, the number of fragments must be on the order of 1/a 2 , 
and so the average number of fragments per channel must be on the order of (in particular, linear 
in) 1/a. As a — > 0, the asymptotics of these estimates become more precise. 
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4. Numbers of Fragments and Gaps 



As mentioned in Section[2J under our probability model, (X(t)) (the process living in the state 
space S) is a Markov process on S. In this section, we obtain analytical results regarding the 
number of fragments and gaps under this process. We begin by formally defining the fragment 
types which were discussed in Section [3] 

Definition 2. For i = 0,1, or 2, a fragment is said to be of type i if it touches exactly i other 
fragments. Ni(t) denotes the number of type i fragments at time t, so that F(t) = No(t) + Ni(t) + 
N2(t) is the total number of fragments. 

Let a(t) denotes the sum of the number of fragments and gaps, 

(1) a{t)=F{t) + G{t). 

The number G(t) of gaps and the numbers of fragment types are related as follows. 
Lemma 1. With probability 1, 

(2) G(t) = N {t) + ±N 1 (t)+I(t) 

For any t > 0, where I(t) = 1, if there is a gap starting at the origin, and otherwise. 

Proof. Each gap, except for boundary gaps starting at or ending at 1, separates two fragments. 
Two gaps surround a type-0 fragment not touching the origin, only one touches a type-1 fragment 
not touching the origin, and none touch a type-2 fragment, so the gaps strictly inside (0, 1) are 
double-counted in 2No(t)+Ni(t). Gaps at the boundaries are counted only once in this expression, 
so if there are gaps touching each boundary, 2 must be added to 2N (t)+Ni (t) to produce a double 
count of all gaps. Then No(t) + Ni(t)/2 + 1 counts the gaps as called for by the lemma. 

With probability 1, a gap always touches the boundary at 1, so the only case left to consider is 
the absence of a gap touching the origin. In this case, there is a type-0 or type-1 fragment touching 
the origin, and so a nonexistent gap has been counted in 2No(t) + N\(t). This over-count cancels 
the under-count of the gap touching 1, and so no correction term is needed, i.e., No(t) + N\(t)/2 
counts all gaps as stated in the lemma. □ 

Definition 3. Let (t k ) denote the sequence of departure times, let Di(t k ) denote the number of 
type i fragments in the channel leaving at time tk, and let A(t k ) denote the number of requests 
admitted to the spectrum at time tk- Finally, define the drift in the total number of fragments and 
gaps: 

(3) Aa(t k ) = f a(t k )-a(t k ^). 

The following lemma is the basis of the stability analysis of a(t) in Section [5j and the 50% rule 
proved later in this section. 

Lemma 2. With probability 1, the departure at t k creates the following change in the total number 
of fragments and gaps: 

(4) Aa(t k ) = A{t k ) - 2D (t k ) - D x {t k ) + J{t k ), k > 1 

with to = 0, and J{t k ) = 1, if a fragment starting at the origin is in the departing channel, and 
otherwise. 

Proof. With probability 1, each new channel allocation covers completely every gap it is allocated, 
except for the last one, which is only partially covered. Thus, with probability 1, each new channel 
allocation changes gaps to fragments, except for the last gap which is changed to a fragment plus 
a gap; this adds one to a(t k -i) for each admission, which accounts for the total of A{t k ) in (j4]). 

Two fragments of the same channel can not be contiguous, so it is correct to add up the changes 
created by departing fragments, with each being treated separately. Suppose first that there is no 
fragment (0, b) against the origin. Then for every type-0 fragment in the departing channel, two 
gaps and a fragment are replaced by a single gap for a net decrease of two, and for every departing 
type-1 fragment, a gap and a fragment are replaced by a single gap for a net reduction of one. 
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This gives the reduction of 2D (t k ) + Di{tk) appearing in If there is a fragment (0, b), it must 
be of type or 1; if it is of type 0, then its departure gives a decrease of one; if it is of type 1, its 
departure has no effect. Each of these contributions is one less than it would be were the fragment 
not touching the origin. There can only be one such fragment, so the correction shown in J(t k ) 
for a fragment (0, b) follows. □ 

We will denote by G~(t k ) the total number of gaps just after the fc-th departure, but before 
new admissions, if any, are made. Note that if we remove A{t k ) from the right-hand side of Q 
and add back the total number of departing fragments at t k , i.e., Do(t k ) + ^i(ife) + -^(ifc), we 
get the number of gaps available to admissions at the k-th departure: 

(5) G-(t fc ) = G(tfc-i) - D Q (t k ) + D 2 (t k ) + J(t k ) 

with J(ifc) = 1, if there is a departing fragment (0, b), and otherwise. 

Knuth's widely known 50% rule appears in a very different context than the model here, so it is 
difficult to anticipate the apparent fact that it also holds for our fragmentation model. However, 
one can argue a similar result assuming only that the fragmentation process has a stationary 
distribution. The result is given below as an expected value of a ratio, rather than a ratio of 
expected values. 

Theorem 1. Assume that the fragmentation process at departure epochs t k has the stationary 
distribution 7r Q; then as a — > 0, 

(G{t k )\ 1 
\R(tk)J 2 



Proof. In the stationary regime, one has, by Lemma [2] 

E^Kifc) -<r(tfc_i)] = 
E^[A(t k ) - 2D (t k ) - Di(t fc ) + J(t fc )] = 0, 

and to balance departure and admission rates, we must have E Ta A(ife) = 1. Thus, we can write 

(6) E ffQ [2£>o(tfc) + £>i(t fc )- J(tfc)] = 1. 

Now for a state x at time t k having r channels and Ni(t k ) type-i (i — 0, 1, 2) fragments, we have 
from Lemma [T] 

E x [2D Q (t k ) + D^tk) - J(t k )} = (2N (t k ) + JVi(t fc ))/r - E x [J(t k )] 

= 2[G(t k ) - I(t k )]/r - E x [J(t k )} - 2G(t k )/r + O (1/r) 

Thus, dropping the 0(l/r) term that tends to with a uniformly in x, the expectation in ([B]) 
proves the theorem. □ 

5. Stability Results 

This section establishes that the average total number of fragments and gaps remains bounded 
and that, for certain distributions of request sizes, crgodicity holds. The analysis leads to the 
two following results. Recall that (t n ) is the sequence of departure times (to = 0) and a(t) = 
F(t) + G(t), defined in {[T|), is the total number of fragments and gaps. 

Theorem 2. There exists some r/ > such that, for any initial state x G S , 



(7) supE^fe^M 

n>l V ' 



< +00. 



Clearly, this implies that for any initial state x € 5, the sequence (E x {a{t n )) 1 n > 0) is bounded. 
With an additional assumption on the distribution of the request size, a stronger stability result 
can be proved. 

Theorem 3. When a > 1/2, the process (X(t)) is positive Harris recurrent; in particular, it has 
a unique stationary distribution. 
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A criterion for finite exponential moments using a Lyapunov function is established next. Then, 
we provide some estimates of the drift of the number of fragments between departures which will 
show us how to construct a Lyapunov function. After constructing this function, we will be in 
position to prove the boundcdncss of exponential moments. 

A Criterion for Finite Exponential Moments. Before stating the main result, some results 
on Markov chains are needed. In the sequel, < st refers to stochastic ordering, i.e., V < s t Z means 
that K(f(V)) < K(f(Z)) for any increasing function /. For reasons that will become clear in 
Lemma the following lemma focuses on admissions at 4 consecutive departure times. Recall 
that (Ui,i > 1) are the sizes of the requests waiting to be allocated bandwidth (after the first one 
U), which are assumed to i.i.d. 

Lemma 3. The random variable A(t{) + ■■■ + A(tj^) is stochastically dominated by a random 
variable Z such that E(e AZ ) < +oo for some A > 0. 

Proof. It is clear that A(ti) < Z + 1 where 

Z = 1 + inf{ra > 1 : C/i + ••• + [/„> 1}. 

The Markov inequality shows that for any z > 0, 

F(Z > z + 1) = F(Ui + ■ ■ ■ + U z < 1) < e (E (e _c/l )) z 

and so E(e r?z ) is finite for r\ > small enough. From this observation, it is not difficult to extend 
the result to A(ti) + ■ ■ ■ + A{t±) instead of just A(t\). □ 

This lemma shows in particular that 

£ = supsupEa;(yl(< 4 )) < +oo. 
i>i xes 

is well-defined; this constant will be used repeatedly throughout the rest of the analysis. The 
proof of the following lemma is standard, and therefore, omitted. 

Lemma 4. Let Z > be a positive, real-valued random variable such that E(e AZ ) < +oo for some 
A > 0, and define c = A~ 2 E(e AZ — 1 — XZ). Then for any < e < A and any real-valued random 
variable V such that V < st Z, we have E (e eV ) < 1 + eE(V) + e 2 c. 

The following result is closely related to result of Hajek [T7]- The proof can be found in 
Appendix [B] 

Theorem 4. Let (Yk) be a discrete-time, continuous state-space Markov chain such that for 
some function f > 0, there exist K, 7 > such that for any initial state y with f{y) > K, 
~ f(Yo)) < ~ 7- Assume that there exists a random variable Z such that for any initial 
state y, Z dominates stochastically the random variable f(Yi) — f{Yo) under W y . Assume finally 
that E(e AZ ) < +00 for some A > 0. Then there exist r\ > and < C < +00 such that for any 
initial state y, 

supE y (e^^A < e^W+C. 

n>l ^ ' 

Theorem |4] will be applied to the Markov chain (X(ti n )) with a function / of the form a K = 
a + nr for some k > suitably chosen. (X(ti n )) is not the most natural choice at first glance, but 
it appears to be needed because of the complexity of the state space. 

It is clear that er(ii) — er(io) < -^(^l) + 1, so that 

<r K (U) - o- K (t ) <{k + l)(A(tx) + ■■■ + A{U)) + 4 

and therefore, by Lemma [3l o~ K {ti) — cr K (£o) is stochastically dominated by some random variable 
Z with an exponential moment. Therefore, one has to establish a negative drift relation for 
f K (^4) — o~ K (to). This is the purpose of the following two subsections. 
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Evolution of the Number of Fragments. Recall that x € S, the initial state of the system, has 
r active channels, and define the total available gap size h = 1— (si+ • • • +s r ). Time referring to 
the initial state x will usually be omitted; e.g., <r(0), A(0), G(0), . . . will be simplified to a, F,G, . . .. 
Recall that Aa(t n ) is defined in ((3j) as a(t n ) — <r(t n _i). 

Lemma 5. Fix < e < 1 and < 77 < 1/2, and let x € S be an initial state such that 
a = G + F > 2K + 1 for some fixed K > 0. 
Then F = N + N t + N 2 > K, and 

1) Ifr = I, then E x (Aa(t 1 )) < £ — K. 

2) Ifr> \ and N Q + Ni > eK , then 

eK 

E x {Aa (h)) <£ + l . 

r 

Assume m i/ie remaining cases that r > 1, define A' ~ K((l — e)/r — e) + , and let i* € {1, . . . , r} 
index a channel Li- in x with the most type-2 fragments. 

3) If Nq + Nx < eK and u> h + Si» , then 



(8) M A <r(t 2 )) < £ 



A' 



r(r - 1) 

4) If iVo + Ni < eK, u < h + s^ and h + Si- < na, then 



(9) E x (Aa(t 3 )) < £• 



r 2 (r — 1) 

5) // Nq + Ni < eK , u < h + Sj» and na < h + Si* , i/ien £/iere exists a 7(77) > suc/i i/iai 

(10) E x (Aa(t 4 )) <C + 2-2l^. 

If follows that there exists a £ > and a function ip(r) > suc/i i/iai /or an?/ 2 wit/i <t > 2A' + 1, 

(11) E x (o-(t 4 )-o-)<l-Ki)(r). 

Proof. As is readily verified, G < F + 1, so 2A + 1 < er = A + G < 2A + 1, and hence F > if as 
claimed. In what follows, we use repeatedly the two following simple facts: 

(12) E x (D (ti) + Di(ti)) = (JV + Ai)/r, 
and by Lemma [1] 

(13) G> A=S> N + N ± > K-l. 

- First case: r = 1. Then, right after the only channel initially present leaves, there is no 
channel allocated bandwidth, and therefore, a(t\) = A(ti). Note that r = 1 is only possible when 
a > 1/2, and in this case the possibility for a channel to be alone is crucial in the proof of the 
Harris recurrence stated in Theorem [3] 

— Second case: r > 1, Nq + N\ > eK. Then the inequality follows from (0}: 

E x (Aa(h)) < £ + 1 - E x (D (h) + D^h)) - $ + 1 - N ° + Nl < £ + l _ 

r r 

In the 3 remaining cases, let N* denote the number of typc-j fragments in any channel i* which 
has the most type-2 fragments. If N0 + N1 < eK, then since F > K, necessarily N2 > (l — e)K and 
-^2 — (1 — e)Kfr. Define the event D* = {channel Li- leaves at ti} and recall that G~ denotes 
the number of gaps right after Li- leaves but before new admissions, if any, are made. It follows 
from ([5]) that G™ > K' in the event D* , since 

G~ = G - N* + N* + J(h) > (-eK + (1 - e)K/r) + = A". 

The remaining analysis tacitly assumes that r > 1, that Nq + N\ < eK and that the channel Li- 
leaves at t\. 
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- Third case: u > h + Sj* . Then -A(ti) = 0, since when Li* leaves it does not provide enough 
additional bandwidth for U . In particular, R{t\) = r — 1 and G(ti) = G~ > K' , and so 

E K (A(T(i 2 )) < £ + 1 - E,(Z? (t 2 ) + D 1 {h); D*). 

The strong Markov property makes it possible to lower-bound this last term. 

E x (D (t 2 ) + £i(i 2 ); D*) = E a (E x(tl) (A)(ii) + Di(*i)); D*) 

= E - Stt V ; ^* > r v *( D ) = 



V J?(ti) ' / ~ r- 1 v ' r(r-l) 
and therefore, E x (Aa(t 2 )) < £ + 2 - K'/(r(r - 1)). 

— Fourth case: u < ft + Sj* < 77a. In this case U is admitted at t\. Thus it makes sense to define 
the event 

E4 = D* n {U leaves at i 2 and f7i > 77a}. 

Then as before 

E x {Aa(t 3 )) < £ + 1 - E x (A,(i 3 ) + Di(t 3 ); £4). 
In the event £4, U is admitted at ii and leaves at t 2 , while U\ stays blocked at and t 2 , so that 
G(i 2 ) = G~ > and R(t 2 ) = r — 1. Hence as in the second case, 

Ex(Do(*3) + -Di(*3);^4) > ^—^(£4) > (1 ~ ^ - 1 

r — 1 r z (r — 1) 

since P^-Ei) = (1 - ij)/r 2 . Thus © holds. 

- Fifth case: u < h + s,;* and i]a < h + s,-«. Again, £/ is admitted at t\. Letting Ui denote the 
sizes of the requests behind U, define the event 

B = {Ui < rja, i = 1, . . . , t and U T +\ > 2r/a} 

with r = inf{n > : U\ H h {/„ > ft + a*. - 77a} and E' 5 = D* n 5 n {[/ leaves at rj 2 }. It is 

readily verified that 1 < r < +00 almost surely. Moreover, one has in E' 5 

< ft* d = ft + s { * - (Ui H 1- U T ) < na < U T+1 . 

This means that at t 2 , exactly r new requests Ui, . . . , U T have been admitted, and U T +± is blocked. 
Moreover, for any i G {1, . . . , r}, one has h* + Ui < 2i]a < f7 T +i, so that if one of the r channels 
allocated to the (fj) leaves, U T +i remains blocked. 

When Li* left, there were G _ > K' gaps; in the remainder of the analysis, we call an initial 
gap a gap present right after L,» left. After Li* left, U and A(t\) — 1 new requests were admitted, 
and then U left and A(t 2 ) new requests were admitted at t 2 . Thus, at t 2 , each initial gap is in 
either of two states: either it is completely filled, or it is still a gap, i.e., it has not been filled 
completely. Let k be the number of initial gaps completely filled at t 2 , and let k' = G~ — k: 
then k + k' = G~ > K'. In each initial gap completely covered at t 2 , there is at least one type-2 
fragment of one of the t new channels. Therefore, N12 + N 2)2 + • • • + AT Tj 2 > k with Ni :2 the 
number of type-2 fragments of the channel corresponding to U. In particular there is a channel 
Lj* , j* G {1, . . . , t} with at least the average k/r of type-2 fragments: Nj* 2 > &/t. Define finally 
the event E5 = E' 5 n {Lj. leaves at t 3 }. Since ft* + L/j* < U T +i, then £/ T +i remains blocked at 
t 3 when E§ occurs, and therefore (note that when j* leaves, some gaps may merge, but not two 
initial gaps), 

G(t 3 ) > Nj* >2 + k' > k/r + k' >(k + k')/r > K' jr. 
Now we proceed as before, to obtain 

E x (A<i{t 4 )) < £ + 1 - E x (D (t 4 ) + D 1 (t i )- E 5 ) 
and, using the Markov property at time t 3 , 

E x (D (U)+D 1 (U);E 5 )=E x (E x(ts) (D (t 1 )+D 1 {t 1 ));E 5 ) 



V R{ta) ' J ~ V r + 
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since R(t^) = r + r — 2 in E§. The same kind of reasoning as before then leads to 



(K' It — 1) + \ K' 
E - ( — - T _2 > E5 ) - -^f(v,h + Si, -a V )-l 



with the function f(rj, •) defined for y > by 

f(T,,y)=E((l + r(y))- 5 ;B(r,,y)) 

with r(y) = inf{n > 1 : Ui + ■ ■ ■ + U n > y} and 

B (v, v) = {Ui < r/a, i = 1, . . . , r(y) and f/ T (i/)+i > ^rja}. 

It is not difficult to show that 7(7?) = inf <j,<i f(rj, y) > which then gives the result. 

It remains to prove (jlip . One only needs to assemble the various bounds, taking into account 
that E x (Aa(ti)) < £ + I for any x E S and i > 0, to arrive at 4 separate bounds on E x {a{t/C) — a). 
For example, using the former bound for the first two terms and the last term of E K (<j(t&) — a) = 
53 1<i<4 Ea;A(T(ti) and then the bound in ([9]) for the third term, we get that one of the 4 bounds, 
which applies when x satisfies the inequalities of the fourth case, is 



(1- V )K' 



Computing the minimum over these bounds with 77 = 1/4 and e — 1/r 2 , one obtains (|11[) after 
setting £ = 4£ + 5 and 

V>« = ^ x ((1 - V) A 7) 
with ip(r) = 1 — 2r~ 2 . This concludes the proof. □ 

We turn now to the case where r is large. In this case, the negative drift comes from the fact 
that, except perhaps at t\, with high probability there is no admission at a departure, since the 
channel that leaves is small with high probability. However, we see in ^ that this is not enough 
for a to decay, for one would need at least one type-0 or type-1 fragment to leave as well, and it 
can be the case that most fragments are of type-2. The second term of the Lyapunov function 
allows us to get around this problem. Since the variation AR(tk) in the number of channels at a 
departure is exactly equal to A(tk) — 1, one readily gets that 

<Tk(U) -o« = (<r(U) -a) + K(A{h) + ■■■ + A{t A ) - 4). 

In particular, if x G S is such that a > 2K + 1 and r < K r , then (from now on, we assume without 
loss of generality that the function ip given by pip in Lemma [5] is decreasing) 

E x (a K {t 4 ) - a K ) <l - K^(K r ) + An{Si - 1) 

whereas if r > K r . , 

E x {a K (t 4 ) - a K ) < 4(C + 1) + K E x (r{U) - r) 
and so we see that we only need to control K x (r(t4) — r) for r large. 

Lemma 6. There exist K r ,^/ r > such that if x € S is such thatr > K r , then E a; (r(<4) — r) < — j r . 

Proof. Since the technical difficulty of the proof of this inequality is similar to that of the above 
proof, we need only give a sketch of it. From S\ H — • s r = 1 — h < 1 one gets #{i : Sj > 7} < 1/7, 
and therefore, P x ( s ii > l/v^ 7 ) < l/v 7 ^ with L il , i\ G {l,...,r}, the channel that leaves at t\. 
Thus, when r is large, with high probability a small channel leaves. 

If h — u is away from a, then the event {u\ > h — u + + ■ ■ ■ + s^ 4 } (with defined similarly) 
has high probability, and in this event A(ti) + ■ ■ ■ + A(t±) < 1. If in contrast h — u is large, then 
with high probability U\ is admitted and with high probability h — u — U\ is away from a; hence 
we can do the same again, and get that, with high probability, A{t\) + ■ ■ ■ + A{t±) < 2. The lemma 
is proved. □ 
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Construction of a Lyapunov Function. For k > 0, one defines, for an initial state x £ S of 

dcf 

the system, a K = a + nr, where r is the number of channels allocated in x and o~ = F + G is the 
sum of the number of fragments and gaps in x. 

Proposition 1. (Lyapunov function inequality) There exist k and K > such that if x € S is 
such that o~ K > K, then E x (cr t .(t4) — er K ) < — 1. 

Proof. Let A r and 7 r be as in Lemma [HI and take n and A as follows: 

4^ + 5 „ 8^ + 2^+2 
k = and A = \- kA,, + 1. 

Assume that a K > K. If r > A" r , then 

E x (cr K (t 4 ) - cr K ) < 4(£ + 1) - n lr = -1. 

Otherwise, r < A r , and since <r K > A , this necessarily gives a > A — fcA" r = 2A" + 1 with 
A = (A — KA" r - l)/2. Thus, 

ExMfc) - O < ? - ^(«r) + 4^ = -1 
and the proposition follows. □ 



Proof of Theorem [2j Theorem 2] and Lemma [3] applied to the Markov chain (X(ti n ), n > 0), 
and the function <r K show that for some r\ > and some constant < C < +oo, 

su P E. T ( e '" T '' ( * 4 " ) ) < e^+C. 

Then the Markov property gives for any i > 

supE,, ( e '" T ' t(t4 "+' ) ) < E x (e^^^ + C < +oo 

from which (|7|) follows readily 



Proof of Theorem [3j In the following discussion, no conceptual argument is missing, only some 
formalism needed to handle the continuous state space S. These details are routine and left to 
the interested reader. In the analysis below, requests are said to be big if their size exceeds 1/2. 

We argue that (X(t)) visits infinitely often a state in which there are no fragmented channels, 
and such that all size distributions remain the same at all visits. This is enough to show Harris 
recurrence; see for instance Asmussen [5]. For this purpose, it is convenient to pick a simple 
regeneration set E C S in which (i) the spectrum is being used by a big request, alone and with 
an unfragmented channel of the form (0, b), and (ii) the request U waiting at the head of the queue 
is also big. Each of these has the conditional request-size distribution given that its size is larger 
than 1/2, i.e., the uniform distribution on (1/2, a), see [22]. 

To verify that E is visited infinitely often, consider the process (R(t), U(t)) with R(t) the number 
of requests allocated a channel at time t and U(t) the size of the request at the head of the queue; 
the process {R{t), U(t)) is simply the process (X(t)) when the data on fragmentation is ignored. 
This process is positive Harris recurrent, as shown in Kipnis and Robert [22] . In particular it visits 
infinitely often states with R{t) = 1 and U(i) > 1/2; one can add U±(t) > 1/2 as well (i.e., the 
first request in line behind the head of the queue request is also big), since this happens with a 
geometric probability. Then when the only channel leaves, the process (X(t)) enters E, since there 
is exactly one channel, it is big, it is necessarily unfragmented and of the form (0,&), and a big 
request is waiting at the head of the queue. Moreover, this argument shows that the time between 
visits to E is integrable, which in turn establishes positive Harris recurrence. This completes the 
proof. 
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Fig. 7. Average number of fragments per channel. 



6. Algorithms 



Although the focus so far has been on measures of fragmentation as a function of a, algorithmic 
issues are also of obvious interest. For example, more uniform patterns of gaps might be an 
advantage. The Linear Scan (LS), discussed in the previous sections, tends to push the gaps 
towards the end of the spectrum, particularly when the spectrum is viewed at random times in 
steady state. Interestingly, our experiments have shown that, for all a < 1/3, the starting position 
of the first gap in the spectrum remains very close to 0.64. 

To uniformize gap locations, an alternative gap scan resembles the Circular Scan (CS) sequences 
of dynamic storage allocation |23j . In our case, CS uses a circular gap list, in which the successor 
to the last gap in [0,1] is the first gap in [0,1]. The scan is still linear, but the starting gap of the 
scan moves as follows: if the last fragment of a channel is placed in gap g, then the residual gap 
of g is the first gap scanned in constructing the next channel. Clearly, although CS will tend to 
uniformize gap sizes as a function of position, boundary effects will persist so long as the spectrum 
itself is not circular, i.e., gaps and fragments are not allowed to overlap the end of the spectrum, 
a restriction that would likely be dictated in practice. 

The average number of fragments per channel is a direct measure of gap-search times, and one 
that we use here. For values of a expected to be of interest in applications, the effects of CS on 
gap-search times are only within a few percent relative to LS, as can be seen in Figure [7J The 
figure also shows the average number of fragments per channel for the Largest-First Scan (LFS) 
algorithm. This algorithm is designed to speed up the process of finding a set of gaps sufficient 
to create a new channel. It selects available gaps in a decreasing order of their sizes and allocates 
them to a request, thereby greedily minimizing the number of gaps needed to fulfill a request. The 
extra mechanism needed for such a search will of course tend to reduce overall performance gains. 
The results in Figure [7] for LFS show a surprisingly large improvement in the average number of 
fragments per channel - as can be seen, a reduction by a factor more than 3 is achieved by LFS 
for even moderately small a. 

The Probability Mass Functions (pmf 's) of the number of fragments per channel are shown in 
Figure HJ Notice that while all probabilities are small under LS and CS, the largest applies to 
the case of no fragmentation at all. The much more peaked distribution for LFS has both much 
smaller mean and variance: The standard deviation under CS and LS is approximately 1.5 to 2.0 
times that under LFS. The same limiting behavior called for by the 50% law holds for CS and 
LFS, which was to be expected, as the arguments supporting the 50% law did not depend on the 
sequence in which gaps were scanned. But an interesting result of our experiments with CS is that 
the 50% approximation to expectation of the ratios is within a couple of percent even when the 
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maximum request size is as much as 0.2 the spectrum size. This is easily seen in Figure |9l The 
convergence rate of LFS to 50% is intermediate between LS and CS. 

7. Normal Approximations 

After the usual scaling (i.e., first centering then normalizing by the standard deviation), the 
scaled version of the number of channels, R, tends in probability to the standard Normal, Af(0, 1), 
as M — > oo, where M = Ll/ a J (it is convenient to express the asymptotics in this section in 
terms of M). This result follows easily from the corresponding heavy-traffic limits in [HE], and 
the ergodicity of (R(t)). 

Interestingly, it was discovered in the experiments that a normal limit law also appears to hold 
for the total number F of fragments as M — > oo (see for example Figure [TU)) . This is not surprising, 
since F is the sum over all requests in system of the numbers Fj,, 1 < i < R, of fragments allocated 
to requests. The requests have a mutual dependence, but one whose effect can be expected to 
weaken for large M. As illustrated in Figure U\ the mean of Fi is proportional to M. Our 
experiment indicate that the standard deviation of Fi is also approximately proportional to M. 
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Fig. 10. Distributions of the total number of fragments and corresponding Nor- 
mal fits. 



Hence, if —¥ denotes convergence in probability as M — s- oo, then 



(14) 



F - f3M 2 
6M 3 / 2 



with p w 1.5 and 9 w 0.9 formalizes the normal-limit law suggested by our experiments. 

Note the two departures from the standard Central Limit Theorem set-up. First, consistent 
with the linear fragmentation observation, individual request fragmentation scales linearly in M, 
so that the total number of fragments scales as M 2 . Restating the limit law in terms of the sums of 
random variables (Fi/M) clearly eliminates this discrepancy. Second, the number of channels (R) 
is random and satisfies the normal limit law discussed above. Thus, the plausibility of (|14|l requires 
an appeal to Central Limit Theorems for random sums (e.g., see |14j . p. 258, for an appropriate 
version). Finally, Figure [10] gives some idea of the convergence of the fits to the normal density; 
as can be seen, fits for a < 0.1 are indeed close to the simulation data. 

We note that the normal approximations shown for the total number of fragments under LS 
were also found to hold under CS and LFS. This is illustrated in Figure ITT1 
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Fig. 11. Distributions of the total number of fragments for a = 0.1 and the 
corresponding Normal fits. 
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Fig. 12. Distributions of the number of gaps at different epochs for a = 0.01. 



The distributions of the number of gaps at different epochs are shown in Figure 1121 The 
distribution at a random time shows a decreasing pmf. What appeared to be yet another normal 
approximation was discovered when looking at the first-admission pmf in Figure 112) this is the 
distribution as seen by the hrst admission immediately after a departure at just those epochs when 
there is at least one admission. The third curve is the pmf of the number, G~ , of gaps as seen 
right after departures, but before determining whether or not the head of the queue fits in the 
total available bandwidth. The fit to the Normal density is also illustrated in the figure. The 
proof of a similar limit law for Renyi's space- filling problem can be found in [9], but extension of 
known techniques once again faces the difficult challenges posed by our more difficult fragmentation 
problem. 



8. Conclusions 

The results of this paper prepare the ground for further research on several fronts. Before listing 
a number of the more important ones, we review what we have learned. Our experiments brought 
out first an unexpected reappearance of a 50% rule relating the expected numbers of gaps and 
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channels in the limit of small request sizes relative to the spectrum size. In our case, we were 
able to prove the limit law. The next result described a linear relationship between the maximum 
request size, a, and the expected number of fragments into which a request was divided at the time 
of allocation. Interestingly, the smaller a was taken, the greater was the resulting fragmentation 
of requests. 

Our stability results established the beginning of a mathematical foundation of fragmentation 
processes. Particularly, we showed that for a > 1/2, the fragmentation process is Harris recurrent. 
For general a, we proved (with considerable effort) that the total number of fragments is bounded 
in expected value. We examined alternative algorithms for sequencing through the available gaps 
and showed that using the LFS algorithm leads to significantly less fragmentation than using the 
LS or CS algorithms. Finally, we exhibited experimentally a limiting, small-a behavior in which, 
with appropriate scaling, distributions tend to Normal. 

A broad direction for further research extends the parameters of our mathematical model. For 
instance, while Uniform distributions are generally the assumption of choice in fragmentation 
models, it would be interesting to see what new effects are created by other distributions of 
request size, e.g., by varying a in the generalized uniform distributions on [0,a], with densities 
x a /a a+1 . The exponential residence-time assumption is likely to yield simplifications to analysis, 
but changes in behavior resulting from other distributions are worth investigating. Moreover, 
instead of a system operating at capacity, in which there is always a request waiting, one could 
adopt an underlying, fully stochastic model of demand, e.g., a Poisson arrival process, as found 
in [22]. 

More realistic, but in all likelihood significantly more difficult models, would relax the inde- 
pendence assumptions. A prime example appropriate for Dynamic Spectrum Access applications 
would be allowing residence times to depend on fragmentation, the greater the fragmentation of 
a request, the longer its residence time. 

The results regarding the performance of the different algorithms imply that the algorithms' 
design should also be considered carefully. Some examples of algorithms that come to mind will 
aim to better fit the fragments into the available gaps. A more challenging objective would be to 
develop algorithms that take into account spectrum sensing capabilities during the gap allocation 
process. 

Finally, another broad and very important avenue of research that introduces more realistic 
models discretizes request sizes and the bandwidth allocation process (as is being done while allo- 
cating OFDM subcarriers). As in other models of fragmentation, the continuous limit represented 
in this paper may conceal important effects, or, conversely, it may introduce effects not present in 
discrete models. We are actively pursuing this avenue of research. 
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Table 1. Averages R, computed from our experiments, vs. the corresponding 
values E(-R), obtained from the numerical computations in [22], for various values 
of a. Except for the case a = 1, the latter results were given to two decimal digits 
only. 



a 


E{R) 


R 


0.05 


39.51 


39.325 


0.1 


19.4 


19.317 


0.15 


12.69 


12.653 


0.2 


9.34 


9.306 


0.25 


7.32 


7.303 


0.3 


5.98 


5.961 


0.35 


5.01 


5.003 


0.4 


4.28 


4.275 


0.45 


3.71 


3.709 


0.5 


3.29 


3.281 



a 


E(R) 


R 


0.55 


2.9 


2.892 


0.6 


2.54 


2.541 


0.65 


2.26 


2.261 


0.7 


2.04 


2.039 


0.75 


1.87 


1.868 


0.8 


1.73 


1.731 


0.85 


1.62 


1.621 


0.9 


1.53 


1.531 


0.95 


1.45 


1.456 


1 


1.392 


1.392 



Appendix A. Tool validation 

In the first test of the simulation tool, we allowed arrivals, discretized the spectrum size, and 
fixed the requested bandwidth, in such a way that the system reduced to a conventional M/M/k 
queue. We then matched simulation results with those obtained from classical formulas. In all 
cases checked, the error compared to the exact results was negligible. 

Keeping with our continuous model operating at capacity with random bandwidth requests, 
one finds that there are very few explicit results for measures of interest, even for the process R{t) 
for the number of channels in the system. Those that do exist can be found in |22j along with the 
elegant formula below, which unfortunately requires calculations that are rarely tractable. 

Theorem 5 (Kipnis and Robert [22)). Let S n denote the sum of n i.i.d. request sizes. Then, the 
maximum throughput (expected number of requests in the system) is 

E(i?) = 



En(l/n)P(S„ < 1 < S n+1 ) E;r 2 (Vn(n - 1))P(S„ > 1) 

A system with requests that are uniformly distributed on [0, 1] is one case that admits of a 
simple formula. Here, the density of S n on [0, 1] is z n ~ 1 / (n — 1)!, and so an easy calculation gives 

E(R) = — !— = 1.392... 

v ' e-2 

Transform methods are used in [22] to evaluate E(R) for other values of a. In particular, com- 
putations are based on the (numerical) inversion of a Fourier-Laplace transform. To evaluate the 
accuracy of our simulator, wc checked our experimental results against the computations for the 
values of a shown in Table [T] As can be seen, the numerics agree well with simulations. 

Appendix B. Proof of Theorem 4 

For any < e < A, 

(15) E y (e^+O) = Ey (E Yn ( e «(/C*)-/W0)) e^l^y^) 

+ E, ( Ey „ (g-C/CttWWO)) e^h {f{Yn]>K} 
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Since by assumption, /(Yi) — /(Yo) under ¥ y is stochastically dominated by Z for every y, one 
gets 

Ey?i ( e *(/(*i)-/C>W) < E ( e eZ) = E ( g AZ) def. ^ 

and therefore 

(16) E y (e y „ ( e -(/(^)-/W»)) e^(^)l {/(y?i) < K} ) < De* K . 

For the second term, we apply Lemma [4] to the random variable /(Yi) — /(Yb) under P^: 

Ey„ (V (/(Yl) - /(yo)) ) < 1 +eEy„(/(Yi) - /(Y )) +ce 2 . 
Thus, on the event {/(Y„) > K}, one gets 

(17) E Yn ( c "(/(*i)-/W,)rt <l- 1£ + ce 2 = f p £ 
and finally 

(18) E, (Ey n ( e *(/C*)-/C«>») e e/(r n ) 1{/(iv)>x} ) < peEy ( e e/(r„)) . 

Gathering the two bounds (|16p and (fT5|) in (|T"5j) finally gives 

Ey (e e/(F " +l) ) < PeE y (e^< V ">) + £e rf 

which leads by induction to 

E v ( e 8 « y »>) < ftfM + De eK . 

Let r\ = (e/(27)) A A. From the definition of p £ in (fT7|) . one can easily check that < p v < 1. 
Therefore, choosing e = r\ gives the result. 
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